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Abstract 

A physical model for the effective thermal conductivity of water frost is proposed for 
application to the full range of frost density. The proposed model builds on the Zehner- 
Schlunder one-dimensional formulation for porous media appropriate for solid-to-fluid 
thermal conductivity ratios less than about 1000. By superposing the effects of mass 
diffusion and eddy convection on stagnant conduction in the fluid, the total effective 
thermal conductivity of frost is shown to be satisfactorily described. It is shown that the 
effects of vapor diffusion and eddy convection on the frost conductivity are of the same 
order. The results also point out that idealization of the frost structure by cylindrical 
inclusions offers a better representation of the effective conductivity of frost as compared 
to spherical inclusions. Satisfactory agreement between the theory and the measurements 
for the effective thermal conductivity of frost is demonstrated for a wide range of frost 
density and frost temperature. 

Key Words: Water frost, effective thermal conductivity, Zehner-Schlunder model, porous 
media, vapor diffusion, eddy convection 
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Nomenclature 


Af s interface area in Zehner-Schlunder model, m 2 

D mass diffusion coefficient, m 2 /s 

h sv latent heat of sublimation, J/kg 

k thermal conductivity, W/m.K 

k e f effective thermal conductivity of frost, W/m.K 

M molecular weight, kg/mol 
p pressure, Pa 

p v partial pressure of water vapor at temperature T s p a 
P v q partial pressure of water vapor at temperature Tq , Pa 

Pr Prandtl number 

R universal gas constant, J/mol.K 
R a gas constant for air, J/kg.K 

R v gas constant for water vapor, J/kg.K 

r radial coordinate, m 
T temperature, K 
T f frost temperature, K 

^ freestream temperature, K 
freestream velocity, m/s 


z axial coordinate, m 


Greek Symbols 

$ solids fraction 

if/ porosity (void fraction) 

v kinematic viscosity, m 2 /s 

X solid-to-fluid thermal conductivity ratio, £2 / ^1 

p f frost density, kg/m 

£ fluid-to-solid thermal conductivity ratio, k\ / ^2 
0) humidity ratio of moist air, kg (water vapor)/kg(dry air) 

AT temperature difference in the enclosure for natural convection, K 

Subscripts 

1 continuous medium (stagnant gas) 

2 solid particle (dispersed) 

00 ambient 

a dry air 
c eddy convection 
d diffusion 
e effective medium 
/ frost 
s frost surface 
t total (ambient) 
v water vapor 
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1. Introduction 

Frost formation represents an important consideration in cryogenics, refrigeration, 
aerospace, meteorology and various process industries. Processes involving simultaneous 
heat transfer and frost formation are occasioned in gas coolers, refrigerators, regenerators, 
freeze-out purification of gases, cryopumping, and the storage of cryogenic liquids. 

In most technical applications frosting is undesirable, except in situations such as 
freeze-drying of foods. Frost offers significant thermal resistance to heat transport in heat 
exchangers, and increases the pressure drop by flow restriction. For example, under certain 
circumstances frost formation can result in 50 to 75% decrease in heat transfer and a 
substantial increase in pressure drop (Emery and Siegel, 1990). Nocturnal frost/icing on 
aircraft wings causes significant aerodynamic penalties in lift and drag (Sahin, 2000), 
leading to a decrease in lift and an increase in drag. Frost/ice formation on cold 
fuel/oxidizer tanks and other components may cause debris concerns when it is shed on 
critical components. A fundamental understanding of the nature of frost formation and 
simultaneous heat and mass transfer during frost growth is thus of great practical interest. 

Frost forms when moist (humid) air comes in contact with a cold surface whose 
temperature is less than the freezing temperature of water (273 K), and is less than the dew 
point temperature, so that water vapor directly passes from a gaseous state to solid state 
(see Burghardt, 1972). The thermodynamic process of frost nucleation is represented with 
the aid of a psychrometric chart (Piucco et al., 2008). 

Among the transport properties of frost (regarded as a porous medium), the effective 
thermal conductivity of frost primarily governs the heat transport through the frost layer 
and its growth rate. Measurements by various investigators suggest that the effective 
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thermal conductivity of frost depends on a number of parameters such as frost density, frost 
temperature, plate surface temperature, and environmental parameters (air temperature T rrj , 
humidity m ca , velocity U x , etc.), the frost density and frost temperature being the most 
important variables. As frost grows on the cold surface and densifies, both the thermal 
conductivity and density increase by an order of magnitude. Several empirical correlations 
have been proposed for the effective thermal conductivity of frost (as a function of frost 
density/or temperature), and form the basis for many theoretical models for predicting frost 
growth. 

An accurate description of the effective frost conductivity is exceedingly difficult. This 
is evident from the nonhomogeneous character of the porous medium, comprised of ice 
particles of various shapes and sizes distributed in moist air, and their interaction involving 
heat and mass transfer. The structure of frost, regarded as a porous medium, consisting of 
different grain (pore) sizes and geometries, thus results in the effective thermal 
conductivity that is far from isotropic. Many models are therefore based on the assumption 
of uniform-sized particles of relatively simple geometry with periodic (regular) structures, 
such as simple cubic (SC), face-centered cubic (FCC) and body-centered cubic (BCC), and 
are generally limited to a narrow range of frost density (solids fraction). For example, the 
maximum values of solid fraction that is possible in a packed bed of uniform spherical 
particles are 0.5236, 0.68 and 0.74 for SC, BCC and FCC types of packing respectively 
(Tien and Vafai 1979). Tien and Vafai (1979) have shown that the SC packing and FCC 
packing respectively yield lower and upper bounds for the effective thermal conductivity of 
a random packed bed. For cylindrical particles, the SC packing yields the limit of 
<fi = 0.7854 . 
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Woodside (1958) proposed a model for the thermal conductivity of snow by 
considering the medium to be composed of a cubic lattice of uniform solid spherical 
particles suspended in a gas, and including the effect of mass diffusion. This model 
considered a unit cube containing one-eighth of a sphere of radius less than or equal to 
unity. Woodside's equation when compared to the thermal conductivity of Pitman and 
Zuckerman (1967) for snow at -88 °C, -27 °C, and -5 °C suggest that it predicts 
considerably steep variation with density for densities in excess of about 400 to 500 kg/m 3 . 
The model is limited to solids fraction of (/> ~ 0.5236 (touching spheres) as a result of the 
unit cube considered. 

Dietenberger (1983) developed a semi-theoretical expression for the frost conductivity 
for the whole range of frost density by postulating a complicated frost structure. A random 
mixture model was proposed with ice cylinders and spheres immersed in moist air at low 
densities, and with ice planes and moist air bubbles at high densities. Mass diffusion is also 
included, and the frost conductivity is shown to depend on both frost density and frost 
temperature. The predicted frost conductivity at high densities undergoes steep variation, 
and considerably (rather unrealistically) exceeds the measurements. 

Le Gall et al. (1987) considered the model of Auracher (1986) for modeling the 
conductivity of frost as a function of both temperature and density . The model is 
applicable for the full range of solids fraction. Apparently unreasonably large values of 
mass diffusion rates (relative to molecular diffusion rates) were postulated to match the 
model with measurements. 

Sahin (2000) developed a model of the effective thermal conductivity of frost during 
the initial crystal growth period, by regarding the frost as formed of cylindrical frost 
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columns (normal to the cold surface) surrounded by moist air, and taking into account mass 
diffusion. Although the model seems to be physically realistic and consistent with observed 
frost structure, its application is limited to early stages of frost formation. 

In most of the theoretical models, the contributions of radiation and convection to the 
effective thermal conductivity of frost have been ignored, and only vapor diffusion 
contribution is considered. While the neglect of thermal radiation and natural convection is 
justified, the role of turbulent eddy convection is open to question. An investigation by 
Biguria and Wenzel (1970) suggest that the turbulent eddy conductivity becomes 
important, and increases in proportion to the freestream velocity. 

In recent times, fractal geometry has been widely applied to predict the thermal 
conductivity of porous media. Yu and Cheng (2002) considered fractal models for the 
effective thermal conductivity of bi-dispersed porous media. The theoretical predictions 
based on fractal models are shown to agree well with experimental data, and the lumped 
parameter models (see Hsu et al. 1994, and Cheng and Hsu, 1999) based on unit cell 
approach. Feng et al. (2007) applied fractal geometry for investigating the effective thermal 
conductivity of unsaturated porous media (3-phases including a solid, a wetting phase such 
as water, and a non-wetting phase such as air) on the basis of self-similarity. 

It is the purpose of this article to propose a theoretical model applicable for the 
effective thermal conductivity of frost over the full range of frost density. The model starts 
from a well-established porous media model devoid of mass diffusion and thermal 
convection effects, and superposes the effect of vapor diffusion and turbulent eddy 
convection on pure conduction in the fluid. The validity of the proposed model will be 
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assessed by comparison with existing correlations and measurements for the frost thermal 
conductivity over an extended range of frost density and frost temperature. 

2. Frost Structure and Mechanism of Frost Formation 

An understanding of the frost structure and the mechanism of frost formation is helpful 
in developing models for the frost conductivity. The frost growth (formation) period can be 
broadly divided into three periods (Cheng and Wu, 2003): 1) a crystal growth period, 2) 
developing frost layer growth period, and 3) fully developed frost layer growth period (Fig. 
la). In the crystal growth period representing early stage of frost formation, one- 
dimensional type ice columns grow out of nucleation sites. In the frost layer growth 
periods, three dimensional growth of ice crystals produce a porous structure of the frost 
layer. The full growth period may be further classified into monotonic and cyclic growth 
periods (Cheng and Wu, 2003). 

Very little is known about the fundamental nature of the nucleation or initial stages of 
frost formation phenomenon; i.e. crystal structure during early stage of frost formation. 

The rapid frost growth rates during the initial stages of frost formation corresponds to the 

/ 

crystal growth period. Crystal structure may be classified as columns, plates and dendrites 
(Sahin, 2000), as shown in Fig. lb. Transitions from a plate growth to prism growth , and 
back to a plate structure and then to prisms again have been observed as the temperature 
decreases. In the crystal growth period of frost formation, the frost is well represented by 
cylindrical frost columns (normal to the surface) surrounded by moist air (Sahin, 2000). 

As frost forms, the frost surface temperature increases above the temperature of the 
cold surface. When the frost surface (air-frost interface) temperature approaches the 
freezing temperature, repeated cycles of melting and freezing occur (Sami and Duong, 
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1989). This results in structural changes in the frost layer that tend to increase both the 
density and thermal conductivity with time without a proportional (corresponding) increase 
in the frost thickness. This phenomenon is more likely to occur in high humidity and/or 
high temperature environments than in low humidity and low temperature environments. In 
the latter case, the frost surface temperature will rise to a steady value below the melting 
point, such that the rate of frost deposition is balanced by the rate of sublimation from the 
frost. Fig. lc displays a general representation of heat and mass transport by diffusion and 
convection in the frost layer. Vapor transport by diffusion is known to be the main 
mechanism by which the temperature dependence of the thermal conductivity is manifested 
(Sturm etal. 1997). 

Frost structure and density exert significant effect on the effective thermal conductivity 
of frost, and depend on frost surface temperature and the surrounding hydrodynamic 
(flow) conditions (Brian et al., 1970). In the low density region when the crystals 
commence to branch out, the convective heat transfer by eddies within the air space is a 
significant contributor to the frost thermal conductivity (Biguria and Wenzel, 1970). For 
example, natural convection frosts formed at low temperatures are usually light and fluffy 
(densities between 10 and 50 kg/m 3 ). Forced convection frosts formed at relatively high 
temperatures (248 - 247 K) have densities ranging from 300 to 500 kg/m 3 , and are 
characterized by relatively high conductivities (Coles, 1954). 

Strictly speaking, the frost conductivity and density are not uniform across the frost 
thickness. In view of the resulting complexity, most of the models proposed for frost 
growth prediction are based on the assumption that the frost density and thermal 
conductivity are uniform throughout the entire frost layer (Le Gall et al., 1997). 
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3. Correlations for Frost Thermal Conductivity 


The first measurements on thermal conductivity of snow were reported by Andrews 
(1886). Since then many investigators reported data on the effective thermal conductivity 
of frost. Unfortunately, there is considerable scatter in the conductivity data arising 
primarily from the differences in the microstructure of snow, with differences in grain 
characteristics (grain size and grain type) and inter-granular bonding (connectivity), leading 
to an order of magnitude range in the conductivity at a given density (Sturm et al., 1997). 
These measurements were presented in the form of correlations of thermal conductivity as 
a function of density and/or temperature. 

A detailed list of various correlations proposed for the effective thermal conductivity of 
frost along with the density limit has been presented in Shin et al. (2003) and Yang and Lee 
(2004). These correlations express the conductivity in terms of frost density and frost 
surface temperature. Most of these correlations are limited to pj < 600 kg/m 3 . 

Van Dusen (1929) proposed the following correlation for the effective thermal 
conductivity of frost: 

k ef = 0.029 + 0.403x10 _3 /? / +0.2367x 10~V} (la) 

where k e j- is in W/(m.K), and pf in kg/m 3 . This correlation was based on the data from six 

other investigators, it does not have a restriction on the range of frost density, but is 
restricted to the average frost temperature to be between 243 K and 273 K. No data on the 
range of air velocities and air humidity is available to the author. This correlation is 
generally known to provide a lower bound for the frost conductivity. Eq. (la) is based on the 
correlation of data of six other investigators. 
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Yonko and Sepsy (1967) suggest the following correlation for the effective frost 
conductivity: 

k e j — 0.024248 + 0.73 lxl 0 -3 /?/ + 0.1 183x1 (lb) 
where k e j is in W/(m.K), and pf in kg/m 3 . This equation is valid for frost density p y < 573 
kg/m 3 , with no lower limit. The data cover a range of humidity co ^ of 0.0075-0.01 5 
(vapor)/kg (dry air), air stream velocity U x of 1 .3-5.3 m/s, and ambient temperature T rxj of 
295 K. 

Ostin and Anderson (1991) proposed the following correlation for the frost thermal 
conductivity: 

k ef = — 8.7 IjcIO - 3 + 4.39xl0~ 4 /> / +1.05x1 0~ 6 /?j. (lc) 

where k e j- is in W/(m.K), and pj in kg/m 3 . This equation is valid for frost density pf in the 
range 50-680 kg/m 3 , T y =253-258 K, and V =3 m/s. 

Other familiar correlations for the frost thermal conductivity as a function of frost density 
were reported by Abels (1893), Jansson (1901), Devaux(1933), Kamei et al. (1950), 
Kandrateva(1954), and Schropp(1935). These correlations have limited range of frost 
density, and are indicated in Yonko and Sepsy (1967), O'Neal and Tree (1985) and Sahin 
(2000). Brian et al. (1970), among others, proposed a correlation for the frost thermal 
conductivity as function of frost density and frost surface temperature, but is limited to pf 

<130 kg/m 3 . 
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The frost density appearing in the correlations for the conductivity is generally correlated 
with frost surface temperature Tj s . For example, Hayashi et al. (1977) proposed the 

following correlation for the frost density: 

p f = 650 exp[0. 227(7^. -273.15)] (2) 

where pj- is in kg/m 3 , and T f s is in degrees Kelvin. Eq. (3) is widely considered for frost 

growth modeling. Eq. (3) is applicable for -25 °C< Tj- s <0 °C, air stream velocities in the 

range of 2-6 m/s, and at an airstream humidity ratio of 0.0075 kg (vapor)/kg (dry air). 

4. Recent Investigations on Frost Thermal Conductivity 

Since the publication of Woodside (1958) model and the measurements of Pitman and 
Zuckerman (1967), many investigators studied the thermal conductivity of snow. Most 
notable among recent experimental work is that of Sturm et al. (1997), which is arguably the 
standard reference for snow thermal conductivity measurements. The correlation proposed by 
Sturm et al. (1997) is expressed by 

k e f = 0. 138-1. OlxlO -3 /?/- + 3 . 233 x 10 " 6 / 9 ^ \56<p f <600 

j j j ( 3 ) 

= 0.023 + 0.234x1 O' 3 /?/ ^ <156 

where p f is in kg/m 3 , and k e p in W/m.K. This correlation is based on very detailed and 

consistent set of data (488 measurements) obtained by the transient method (short test 
duration) with cylindrical samples (of various types of snow) placed in an environmental 
chamber (to control the sample temperature). 

Recently, the thermal conductivity of snow has been modeled by Kaempfer et al. (2005) 
as a complex porous media with the aid of three-dimensional geometries of snow structure. 
Kaempfer et al. (2005) also compared their model thermal conductivity with a single point 
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measurement at a density of 268 kg/m 3 (</> = 0.29 ). More complex snow geometries were 
incorporated in the models of Aron and Colbeck (1995, 1998). 

In this manner, much work has been done to improve the models of thermal conductivity 
through microstructure, well beyond assuming a simplified geometry of spheres. 

5. Proposed Model 

In the proposed theory, we regard the frost as a porous medium. We consider as a starting 
point the model of Zehner and Schlunder(1970) established for packed beds, which is known 
to be valid for the full range of void fractions and for solid-to-fluid conductivity ratio less 
than 1000 (Cheng and Hsu, 1999; see also Kaviany 1995). As the ice to air thermal 
conductivity ratio is about 100, the choice of Zehner-Schlunder model is justified for our 
current purposes. Vapor mass diffusion and forced convection effects are then taken into 
account to arrive at the effective frost thermal conductivity. 

The following physical assumptions are made: 

(a) Heat transfer in the frost layer is primarily one-dimensional. 

(b) Frost density and thermal conductivity are independent of frost thickness. 

(c) Thermal radiation within the frost layer is negligible This is so considering that the frost 
temperatures are low. Dietenberger (1983) has shown that for ice crystals the radiation 
effective conductivity is negligible for the size of ice crystals and the temperature levels of 
the frost layer (an order of magnitude less than the thermal conductivity of air). 

(d) The gas spaces in the frost layer are small enough that heat transfer by natural 
convection within the frost layer is negligible. The effect of natural convection in porous 
frost layers is manifest on account of convective currents (Benard convection represented by 
a pattern of hexagonal cells) set up in enclosed spaces under certain conditions of unstable 
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configuration. Natural convection becomes significant only when the gas-filled spacing is of 
the order of 1 cm (Woodside, 1958). 

(e) The total gas pressure is constant within the frost layer, and is equal to the atmospheric 
pressure (Le Gall et al., 1997). 

(f) Soret and Dufour effects are negligible (Le Gall et al., 1997). 

(g) The local humidity of air in the frost layer corresponds to saturation humidity at that 
temperature (thermodynamic equilibrium). 

A description of the various components of the effective thermal conductivity is provided 
below. 

5.1. Packed Bed Stagnant Thermal Conductivity 

Many theoretical models for the effective thermal conductivity of porous packed beds 
were developed since the classical work of Maxwell (1873) and of Lord Rayleigh (1892), 
which were limited to relatively small solids fraction. 

Zehner and Schlunder (1970) arrived at an analytical expression for the effective stagnant 
thermal conductivity of a packed bed . They considered for the unit cell one-eighth of a 
cylinder (inner cylinder of unit radius and outer cylinder of radius R ), see Fig 2. Fluid is 
filled between the inner and outer cylinders while the inner cylinder consists of both the solid 
phase (shaded area) and the fluid phase with the interface area Aj s described by (Cheng and 


Hsu’ 1999) 


r 2 +- 




= 1 


(4a) 


where r and z refer to the cylindrical coordinates (Fig. 2), and B is the shape factor 
characterizing the geometrical effect of the solid particle. For B — » 0 , the boundary becomes 
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the z axis with no solid volume; for 5 = 1 the solid becomes a sphere, and for 5 — » cothe 
solid occupies the entire inner cylinder. In the cylindrical geometry, the cylinders are 
oriented vertically in the direction of heat flux. 

In their analysis, Zehner and Schlunder(1970) assume that one-dimensional heat 
conduction in the unit cell (Fig. 2) follows two parallel paths: the first path is through the 
fluid in the outer concentric cylinder (of inner radius of unity and outer radius R ), and the 
second path is through the inner cylinder (of unit radius), which is comprised of both the 
solid and the fluid phases. Thus the effective thermal conductivity of the unit cell is 
expressed by 


f 


1 






+ 


l 


ei 


R 2 k x 


(4b) 


The quantity k ei is the equivalent thermal conductivity of the inner cylinder, and depends on 


the shape of the solid-fluid interface A j s given by (Eq. 4a). Based on an analogy with mass 


transfer experiments, a relation for the unknown R in terms of the porosity is considered as 

MR 2 = yjl-ys 


(4c) 


For a more detailed derivation of the effective thermal conductivity, the reader is referred to 
the work of Hsu et al. (1994). 

The effective thermal conductivity of a porous packed bed is finally expressed by 
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V 


- = l-JT^y/ + 
k\ 1 L(l-£5) 2 


(1 


In 




K&J 


5 + 1 


V ^ J 


5-1 
1 -4B 


(5a) 


where the shape factor 5 is related to the porosity y/ by the relation 


5(3 - 45 + 5 2 + 21n5) ~[ 2 


(5 - 1) 3 


= 1-^ 


(5b) 
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and 


% = k\ /&2 =1 IX (5c) 

In the above equations, k\ and ko denote the thermal conductivity of the stagnant fluid and 
the solid respectively. The porosity y/ appearing in Eqs. (5a) and (5b) is defined by 

¥ = {P2-P)I{P2~P\) ( 5d ) 

where p stands for the density of the packed bed (frost in our case), p\ is the density of the 
fluid (air-vapor in our case), and p 2 the density of the solid (ice in our case). It is related to 
the solids fraction <j> by the relation 


y/ = \-(j) 


(5e) 


It was suggested by Zehner and Schlunder(1970) that Eq. (5b) can be approximated by 


B = C 




(5f) 


^ V ) 

with m = 1 0/9 . The constant C depends on the shape of the particle. Zehner and Schlunder 
(1970) suggest that C = 1.25, 1.40 and 2.5 for spheres, broken (irregular particles) and 
cylinders respectively. 

Cheng and Hsu (1999) point out that the Zehner-Schlunder model (1970) agrees well with 
the experimental data of Nozad et al. (1985) at 5 / = 0.4, involving spherical inclusions for 


&2 / k\ < 10 . For kj / k\ > 10 , the Zehner-Schlunder model underpredicts the effective 

stagnant thermal conductivity substantially. Hsu et al. (1994) postulated that for / k\ >10 
solid-solid conduction through the finite contact area between the spherical particles (arising 
from flattening of spheres due to external loads or the weight of the bed itself) becomes the 
main mechanism of heat transfer, and that the reason for this underprediction of data for 
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kj !k\ > 10 is the assumption of point contacts between spheres. For values of kj / k\ < 10 , 
heat conduction through the finite contact area is of lesser importance, and heat transport is 
primarily controlled by the particle-fluid-particle conduction. A detailed comparison of 
Zehner-Schlunder model with data for k^ / k\ ranging from 8 to 1200 and for various values 
of solids fraction was recently presented by the author (2010). 

5.2. Effect of Mass Diffusion 

The effective thermal conductivity of air in the presence of vapor diffusion can be 
expressed as 

k \d = k \ +k d ( 6 ) 

where k d is the contribution of mass diffusion to the effective thermal conductivity of air. 

The contribution of thermal radiation to the effective conductivity of air is neglected relative 
to the effects of diffusion (Dietenberger 1983). 

The mass diffusion contribution kj is determined from (Sahin 2000) 

, . _ dco 

k d = h sv D Pa ~=r ( 7 ) 

dl j 


where co denotes the local humidity ratio of moist air in the frost layer (assuming 
equilibrium), h sv the latent heat of sublimation, D the mass diffusivity of water vapor. We 
can rewrite Eq. (7) in the form 


k d 


k svDpa 


dco dpg 
dp g dTf 


( 8 ) 


where 


co = 0.622 p v / p a 


( 9 ) 
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Pa ~ Pt ~Pv 


( 10 ) 


In the above equations, p t is the ambient (total) pressure (taken as 1 atm here), p v is the 
partial pressure of water vapor, and p a the partial pressure of dry air. The mass diffusivity is 
determined from (Pruppachar and Eklett , 1978) as 

D = 2.11xlO“ 5 (r / /r 0 ) 1 - 94 m 2 /s (11) 

where T j- is the frost temperature in degrees Kelvin, and 7 q the reference temperature 
(273.15 K). The partial pressure p v at temperature Tf is given by the Clasius-Clapeyron 
equation 


Pv=Pv oexp— — (12) 
p v yjQ if j- 

where p v Q stands for the partial pressure of water vapor at reference temperature Tq , and R v 
the gas constant for water vapor. Noting that 

Pa~P,faTf) 03 ) 

we express Eq. (8) in the form 


, , 2 D(M a /M v h sv 

kd = K T ~T Pv 0 ex P V - 

RvT 3 Ry 


where M a , M j denote the molecular weights of air and water vapor respectively. 

5.3. Effect of Eddy Convection 

Most of the data on frost thermal conductivity measurements were obtained for flow of 
moist air under forced convection, while the remainder belong to natural convection. At 
relatively higher densities the frost becomes a mesh of dendrites with enclosed air pockets 
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(Biguria and Wenzel 1970). It is therefore plausible that a component of this air flow is 
directed transverse to the frost layer composed of surface roughness elements, thus randomly 
penetrating the voids and enhancing the energy transport by turbulent eddy convection. 

In the presence of mass diffusion and eddy convection, the effective thermal conductivity 
of air may be expressed as 

k \dc =k\+k d +k c (15) 

where k c denotes the contribution of eddy convection to the effective conductivity of air. 

By a comparison of their experimental data for frost conductivity and analysis, Biguria 
and Wenzel (1970) have arrived at a correlation for the turbulent (eddy) conductivity of air as 
a function of freestream velocity U 0 0 as follows: 

k c = 0.00568 1/* W/(m.K) (16) 

where U ^ is expressed in m/s. According to the data of various investigators compiled by 
Yang and Lee (2004), forced convection measurements were obtained for freestream 
velocities ranging from lto 12 m/s. The data compiled by O'Neal and Tree (1985) for 
cylinders, flat plates, annuli, and parallel plates suggest that the freestream velocities ranged 
from 0.2 to 17 m/s. 

Equation (16) indicates that k c I k\ =1.2, 2.4 and 3.5 for U x =5, 10 and 15 m/s 
respectively, considering k\ ~ 0.024 W/(m.K). This consideration represents the significance 
of eddy convection in enhancing the effective thermal conductivity of frost by increasing the 
effective stagnant conductivity of air. This observation explains the considerable scatter in 
many of the correlations proposed for the effective thermal conductivity of frost as a function 
of frost density alone. 
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From a more fundamental consideration, the eddy convection in porous media is 
intimately related to the phenomenon of thermal dispersion (see Cheng and Vortmeyer, 
1988; Hsu and Cheng, 1990). Thermal dispersion is manifested on account of microscopic 
fluid velocity and temperature deviations from the corresponding average values, and results 
in enhanced heat transport. Thus fluid flow in the pore channels of porous media my be 
regarded as pseudo-turbulence. Many correlations were proposed for transverse and 
longitudinal thermal dispersion conductivity, and substantial variation prevails between the 
various correlations (Du et al., 2003). Experiments on forced convection through porous 
media suggest that the magnitude of transverse thermal dispersion conductivity is often much 
larger than that of molecular diffusion, and that the axial thermal dispersion conductivity is 
negligible relative to the bulk convection (Yu and Li 2004). Measurements also reveal that 
thermal dispersion increases as the porosity decreases (Du et al., 2003). 

Hsu and Cheng (1990) proposed correlations for the transverse thermal dispersion 
conductivity for high and low Reynolds number (Re y ) flows in porous media as: 


where 

and 



0 -^) 

¥ 


Re f » 10 


Re y 


= D 2 ^—j^-Pe 2 f Re y « 10 
y/ 2 J 

Pey = Rey Pry 
= u m d p lcif, P r / = v f ! a f 


(17a) 


(17b) 

(17c) 


with the subscript denoting the fluid phase. The quantities Pry ,Pep denote respectively the 
Prandtl number and Peclet number of the fluid, u m the mean axial velocity of the fluid, and 
d p the particle diameter. The constants D \ , D 2 are evaluated by comparison of the model 
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with measurements. The effect of tortuous flow paths on the constant Z)j for high Reynolds 
number flow was reported by Yu and Li (2004) with the aid of a fractal model. Transverse 
thermal dispersion and wall channeling effects (existence of peak velocity in a region close to 
the external boundary) were treated by Cheng and Vortmeyer (1988). Convection in variable 
porosity media, including channeling effects, were considered by Vafai (1984). 

In the present context of frost growth, the external flow is on the top of the frost layer 
instead of streaming through the interior of the frost layer in the axial direction transverse to 
the frost thickness, and only a portion of this external flow is directed into the interior of the 
frost layer in longitudinal and transverse directions, leading to low Reynolds number flow 
regime. Thus existing information on the thermal dispersion effects in the traditional packed 
beds cannot be directly applied to the circumstances of frost formation. Much work remains 
to be undertaken before quantitative effects of thermal dispersion in frost layers can be 
addressed with regard to the effective thermal conductivity of frost. A detailed account of 
thermal dispersion effects in frost layers is beyond the scope of the present work. 

5.4. Thermophysical Properties 

In the present analysis, the properties are taken as follows. The thermal conductivity of 
air is calculated as (Rohsenow and Hartnett’ 1973; Dietenberger, 1983) 


k\ ~ ^air = 2.646x10 


-13 


4/2 


[l + (245/r)l0 


12 IT 


w 

m.K 


(18) 


where T is in degrees Kelvin. The thermal conductivity of ice is obtained from the relation 
(Dillard and Timmerhaus, 1969; also Dietenberger, 1983) 

k 2 =k ice =630/T W/( m. K) (19) 

where T is the absolute temperature (degrees K). 
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The latent heat of sublimation of ice is determined from (Mago et al., 2004) 

hfr ={-0.1083[l.8(r-273.16) + 32]+2833} kJ/kg (20) 

where T is in degrees Kelvin. 

6. Results and Comparison 
6.1 Comparison with Correlations 

In the evolving frost layer, the actual structure of the ice particles is very complex, and 
far from being that of spheres. In an effort to understand the sensitivity of the effective 
conductivity of porous media on particle shape, calculations with Eq. (5) have been made for 
spheres, cylinders and irregular particles with k^lk ^- 100 in the absence of mass diffusion 
and eddy convection, and the predicted effective thermal conductivity is compared with the 
correlations of frost thermal conductivity to density due to Yonko and Sepsy (1967) and of 
Van Dusen (1929). The results suggest that cylinders lead to the highest effective 
conductivity with the irregular particles producing slightly higher values than the spheres, 
and that the cylinder results are closer to the experimental data (correlations) of frost relative 
to those of spheres and irregular-shaped particles. In view of these considerations, we 
hypothesize here that the frost particles are better represented by cylinders, and thus we 
consider only cylinders in our subsequent results. 

In the case of cylindrical geometry referred to in the Zehner-Schlunder model, the 
cylinders are oriented vertically (normal to the plate surface) in the direction of heat flow. 
This choice of cylinders points to a very real effect in the formation of frost, which can 
develop anisotropically, with a preferred direction of growth. This consideration might 
explain why the cylindrical geometry is better suited over spheres in describing the frost 
thermal conductivity. 
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The role of mass diffusion on the frost thermal conductivity is illustrated in Fig. 3, with 
the aid of Eqs. (5), (6), and (14), with Arj in Eq. (5a) replaced by . The mass diffusion 

effect is shown for three different frost temperatures of 250 K, 260 K, and 273 K (the latter 
representing the maximum possible value). The conductivity contribution due to mass 
diffusion increases with an increase in temperature. There is seen a crossover of the curves at 
large value of <p . At Tj ■ = 273 K, the mass diffusion contribution improves the agreement 

with data; but even with mass diffusion, there still exists considerable underprediction of the 
data. Thus additional mechanisms for enhancing the frost conductivity seem to exist. 

Figure 4 presents the predictions of k e p with mass diffusion and eddy convection on the 

basis of Eqs. (5) and (14-15), with k\ in Eq. (5a) replaced by k\j c . The diffusion 
contribution is taken at 7y = 273 K, and an eddy convection factor k c / k\ = 1 is considered. 

The increase in the conductivity due to eddy convection is of the same order as that due to 
mass diffusion. In the Yonko and Sepsy 's data, the air velocities range of 1 .3 to 5.3 m/s 
roughly correspond to k c / k\ of about 0.3 to 1 .25 (an average of about 0.8) in accordance 

with Eq. (16). With the inclusion of eddy convection effect, the agreement between the 
theory and the data has improved. In particular, the theory accounting for combined diffusion 
and convection matches closely with Van Dusen's correlation up to about ^ = 0.4 , but still 
significantly underpredicts the conductivity expressed by the Yonko-Sepsi correlation. A 
frost temperature of 273 K is slightly unrealistic in terms of experimental considerations, and 
so are the conditions of natural convection in the experiments. 

The effect of eddy convection levels on the effective thermal conductivity of frost, as 
computed from Eqs. (5) and (14-15), is demonstrated in Fig. 5. The comparisons are shown 
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for k c I k\= 1,2 and 3, and at T = ITS K. It is observed that with increasing value of the eddy 
convection, the effective frost conductivity increases, as is to be expected. Recall that the 
Yonko-Sepsy correlation constitute data for which k c / k\ ranges from 0.3 to 1.25. For 
(j> <0.6, better agreement is achieved between the prediction and the Yonko-Sepsy 
correlation with a choice of k c / ky - 3 . In the case of Van Dusen correlation, the condition 
of k c lk\= 1 affords better accord between the theory and the correlation up to (j) = 0.4 . It 
appears that a mean value of k c / k\ = 2 may be regarded as reasonable considering the lull 
range of frost density. It appears that the value of k c / k\ = 3 represent the high end of free 
stream velocity values, and are not as realistic as k c / k\= 1 or 2. 

A comparison of the present theory, according to Eqs. (5) and (14-15), with the 
correlation of Sturm et al. (1997) is depicted in Fig. 6. The data are obtained for mean 
temperature of about 258 K. cover various types of snow, having different microstructure, and 
as a result show considerable scatter. The results from Van Deusen (1929) correlation and 
from Ostin and Anderson (1991) correlation are also shown for comparison purposes. For 
solid fractions ^between 0.2 and 0.5, the theory is closer to Van Deusen and Ostin and 
Anderson correlations, whereas for (ft > 0.5 the theory agrees better with the correlation of 
Sturm et al. (1997) and of Ostin and Anderson correlation, which extends up to the farthest 
value of ^ « 0.74 . 

6.2 Comparison with Measurements 

Figure 7a illustrates a comparison of the present theory, given by Eqs. (5) and (14-15), 
with the data for the frost conductivity obtained by Pitman and Zuckerman (1967) and 
Kaempfer et al. (2005). Eq. (18) is used to obtain the thermal conductivity of air. Both the 
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measurements correspond to an average temperature of 246 K. With regard to the data of 
Pitman and Zuckerman (1967), the snow used in the experiments was composed of vapor 
grown crystals having length of 3 to 5 mm. The estimated uncertainty in the measured 
conductivity is about 10 percent. In the measurements of Kaempfer et al. (2005) the snow 
sample is 2 cm high and 4.8 cm in diameter, with the data taken at a single density of 268 
kg/m 3 . The measured conductivity is 0.18 ±0.02 M/m.K. The snow microstructure was 
imaged by X-ray micro-tomography. The predicted snow conductivity from their 
microstructural numerical model (based on some assumed structural parameters, and 
accounting for conduction mechanism only) was found to be 0.15 ± 0.01 W/m.K. 

The comparisons suggest that the present theory, accounting for mass diffusion and an 
eddy convection factor k c / k\ = 1 , matches the data satisfactorily over a broad range of solid 

fraction up 0 = 0.45 (density of 450 kg/m 3 ). Only at 0 = 0.54 the theory with k c / k\ =2 

offers a better agreement with the data. The relatively close agreement between the two sets 
of data at 0 = 0.29 (density of 268 kg/m 3 ) should serve as a testament to the reliability of the 
data over the particular range of frost density. These comparisons suggest the importance of 
mass diffusion and eddy convection phenomena in governing the effective thermal 
conductivity of frost. 

Figure 7b shows a comparison of the predicted frost conductivity, provided by Eqs. (5) 
and (14-15), with the data of Pitman and Zuckerman (1967) for temperatures of 185 K, 246 
K and 268 K. These data are the most widely cited data, and were obtained with an artificial 
snow (freezer frost) with a steady state measuring method. Note that in the ordinate the 
dependence of air thermal conductivity with temperature &i(r)is taken into account. The 
results suggest that the data for thermal conductivity at different temperatures (covering a 
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wide temperature range) when presented in dimensionless form nearly collapse up to <j> of 
about 0.45, and in this range excellent agreement between the theory and the data is noticed. 
For ^>0.5, the measured conductivity at the highest temperature of 268 K considerably 
exceeds the predictions. According to Sturm et al. (1997), the reliability of the measurements 
of Pitman and Zuckerman (1967) is questionable. 

The comparisons presented in Figs. 7a and 7b, suggesting that the model underpredicts 
the data for the effective frost thermal conductivity for ^exceeding about 0.5, reveal that 
thermal dispersion effect likely increases with solids fraction in this region. This is 
qualitatively consistent with the measurements on packed beds suggesting that thermal 
dispersion effect increases as the porosity decreases (Du et al., 2003). 

A comparison of the model predictions from Eqs. (5) and (14-15) with the experimental 
data of Brian et al. (1970) for the average thermal conductivity at Tj = 175K is shown in 

Fig. 8. Thermal conductivity of air is obtained from Eq. (18). The data are limited to 
relatively low density frost having density ranging from 50 to 130 kg/m 3 , The conductivity 
data are obtained from transient measurements. The average thermal conductivity of the frost 
layer was determined from the measured values of the heat flux, the average temperature 
difference across the frost, and the frost thickness (Fourier's law). The scatter in the data is 
believed to be due to structural differences in the frosts of the same density, and due to 
differences in the freestream velocity and humidity. The present theory based on mass 
diffusion and an eddy convection factor of one agrees favorably with the measurements. 
These comparisons, although restricted to a relatively lower range of frost density, appear to 
be consistent with those presented in Fig. 7, and serve to support the validity and usefulness 
of the proposed model. 
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7. Discussion 


The present results suggest that cylindrical inclusions may better represent the frost inner 
structure. As the heat transfer in frost is anisotropic, the model presented here is not expected 
to be applicable to the prediction of the effective thermal conductivity of frost in the direction 
perpendicular to the main axes of the cylindrical inclusions. However in most of the one- 
dimensional formulations of heat and mass transfer during frost growth, it is the effective 
thermal conductivity of frost in the direction of heat flow that is required. A detailed account 
of two-dimensional thermal conduction in frost layers is beyond the scope of the present 
investigation. 

8. Conclusions 

The Zehner-Schlunder model for the effective thermal conductivity of porous media has 
been extended to describe the effective thermal conductivity of frost. By superposing the 
physical effects of mass diffusion and eddy convection on the molecular thermal conductivity 
of air, as expressed by Eqs. (14) and (15), the proposed model yields estimates for the 
effective thermal conductivity of frost, given by Eq. (5), as a function of frost density and 
temperature over the full range of frost density. Comparisons of the model predictions with 
the well-known correlations and measurements for frost conductivity suggest that the frost 
may be better represented by an assemblage of cylinders rather than spheres. The proposed 
model reveals that both mass diffusion and eddy convection effects have significant effect 
on the frost thermal conductivity, and are of the same order. 
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Figure Captions 


Figure 1(a). Representation of frost structure during various stages of frost growth. 

Figure 1(b). Representation of ice crystal structure (Sahin, 2000). 

Figure 1(c). Representation of heat and mass transport by diffusion and convection 
in a frost layer. 

Figure 2. Schematic of the unit cell in Zehner-Schlunder model (1970). 

Figure 3. Effect of mass diffusion on the effective thermal conductivity of frost at various 
frost temperatures. 

Figure 4. Effect of mass diffusion and eddy convection on the effective thermal conductivity 
of frost at 273 K. 

Figure 5. Dependence of thermal conductivity of frost on mass diffusion and eddy convection 
for various eddy convection coefficients. 

Figure 6. Comparisons of predicted thermal conductivity of frost with the correlations of 
Sturm et al. (1997), Van Deusen (1929), and Ostin and Anderson (1991). 

Figure 7a. Comparisons of predicted thermal conductivity of frost with the measurements of 
Pitman and Zuckerman (1967) and Kaempfer et al. (2005) at 246 K. 

Figure 7b. Comparisons of predicted thermal conductivity of frost with the measurements of 
Pitman and Zuckerman (1967) and Kaempfer et al. (2005) at various temperatures. 

Figure 8. Comparisons of predicted thermal conductivity of frost with the measurements of 
Brian et al (1970). 


33 


-» time 
nuclcation 


conden- 

sation 


Y 


fully-developed 

growth 


l-d growth 
" frost column 


n r lair 

h ) i rru'i y 7 ^ ; ■ 

supercooled frost nucleus cooled surface 
droplet 




Figure 1. (a) Representation of frost structure during various stages of frost growth, (b) 
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Figure 2. Schematic of the unit cell in Zehner-Schlunder model (1970). 
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Figure 3. Effect of mass diffusion on the effective thermal conductivity of frost at various 
temperatures. 



<t> 

Figure 4. Effect of mass diffusion and eddy convection on the effective thermal conductivity 
of frost at 273 K. 
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Figure 5. Dependence of thermal conductivity of frost on mass diffusion and eddy convection 
for various eddy convection coefficients. 
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Figure 6. Comparisons of predicted thermal conductivity of frost with the correlations 
of Sturm et al. (1997), Van Deusen (1929), and Ostin and Anderson (1991). 
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Figure 7a. Comparisons of predicted thermal conductivity of frost with the measurements of 
Pitman and Zuckerman (1967) and Kaempfer et al. (2005) at 246 K 



♦ 

Figure 7b. Comparisons of predicted thermal conductivity of frost with the measurements of 
Pitman and Zuckerman (1967) and Kaempfer et al. (2005) at various temperatures. 
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Figure 8. Comparisons of predicted thermal conductivity of frost with the measurements of 
Brian et al (1970) at various temperatures. 
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